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Abstract. Unexpected relations are found between the Toda lattice, the relativistic Toda 
lattice and the Bruschi-Ragnisco lattice, as well as between their integrable discretizations. 



1 Introduction 



We continue in this paper playing some curious game named "reparametrizations" with 
the well-known integrable lattice equations of the classical mechanics. This game was 
started in the papers 0, where several "new" integrable lattices were introduced 
together with their integrable discretizations. 

In the present paper we shall demonstrate that the equations of motion of one of 
simplest flow of the relativistic Toda hierarchy may be brought into the form: 

Xk = Xk [Xk+i - 2xk + Xk-i) + . (1.1) 

\Xk+l-Xk Xk-Xk-l) 

This form of equations of motion is very remarkable, since the right-hand side may 
be seen as a sum of the right-hand sides of another two well-known integrable lattices. 
The flrst one of them is the usual Toda lattice, the equations of motion for which can be 
put into the form: 

Xk = Xk{xk+i-2xk + Xk~i) ■ (1.2) 
The second one is exactly the so called Bruschi-Ragnisco lattice 0, 0: 

Xk+lXk XkXk—1 „N 

Xk = . (1.3) 

Xk+l ~ Xk Xk — Xk-l 

Even more bizarre becomes the situation when looking at the discrete time coun- 
terparts of these systems. Namely, we shall demonstrate that the discretization of the 
relativistic Toda flow derived in in the parametrization corresponding to ( |1 . 1| ) takes 
the form 



. ^ r ~ \ (l + h{xk+i -Xk] 

Xk-Xk [Xk+l-Xk) [Xk-Xk-lj ^ 



Xk Xk 



(Xk+l-Xk) (Xk-Xk-l) (l + h{Xk - Xk^l)) 



;i.4) 



The right-hand side of this equation may be seen as a product of the right-hand sides 
of the following two equations: 

^ _ 1 + h{xk+i-Xk) 

^-^^ = T-TJ——v (1-5) 

Xk- Xk 1 + h[xk - Xk-i) 



and 



Xk Xk 



{Xk+l - Xk) [Xk - Xk-1) 



Xk — Xk {Xk+1 —Xk) {Xk — Xk^i) 



;i.6) 
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Here ( |1.5|) turns out to be exactly the same discretization of the Toda lattice as 
derived in |Q , ||^ and discussed further in , but expressed in the parametriztion leading 
to (|1.2|) . And ( |1.6|) is an integrable discretization of the Bruschi-Ragnisco lattice which, 



to the author's knowledge, has not yet appeared in the literature, and can be derived on 
the basis of a general recipe worked out in |^, |Q and the subsequent publications of the 
author. 

All the systems above (continuous and discrete time ones) may be considered either 
on an infinite lattice {k E Z), or on a finite one {1 < k < N). In the last case one of the 
two types of boundary conditions may be imposed: open-end {xq = oo, Xtv+i = — oo) or 
periodic (xq = x^, x^+i = Xi). We shall be concerned only with the finite lattices here, 
consideration of the infinite ones being to a large extent similar. 

2 Simplest flow of the Toda hierarchy 

We consider in this section the simplest flow of the Toda hierarchy. All the results here are 
not new but collected in the form convenient for our present purposes. For the relevant 
references see [0. 

The simplest flow of the Toda hierarchy is: 

hk = akipk+i - bk) , bk = ak-ak-i. (2.1) 

It may be considered either under open-end boundary conditions (oq = ctAr = 0), or under 
periodic ones (all the subscripts are taken (mod A^), so that oq = un, bN+i = bi). 
The flow ( p.l|) is Hamiltonian with respect to the following Poisson bracket: 



{ctfc, bk} = -{ak, bk+i} = ak (2.2) 
(only the non-vanishing brackets are written down), with the Hamiltonian function 

1 N N 

HT = -Y.bl + T.^^- (2-3) 

^ k=i k=i 

An integrable discretization of the flow ( p.l|) is given by the difference equations 
ak = ak——, bk = bk + h\ — -- — , (2.4) 

Pk \Pk Pk-l J 

where (3k = A (a, b) is defined as a unique set of functions satisfying the recurrent relation 

Pk = l + hbk- (2.5) 
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together with an asymptotic relation 



Pk = l + hbk + 0{h^). (2.6) 



In the open-end case, due to = 0; obtain from ( |2.5| ) the following finite continued 
fractions expressions for 

Pi = l + hbr, P2 = i + hb2 • 



1 + hbi 



(3n = 1 + hbN 72 

1 + /lOjv„i 



1 + hbN-2 



1 + hbi 

In the periodic case (|2.5| ), ( p.6|) uniquely define PkS as A^-periodic infinite continued 
fractions. It can be proved that for h small enough these continued fractions converge 
and their values satisfy ( p.6|) . 

It can be proved ^ that the map ( [2.4] ) is Poisson with respect to the bracket ( ^.21 ). 

Let us recall also the Lax representations of the flow and of the map ( |2.4| ). They 
are given in terms of the N x N Lax matrix T depending on the phase space coordinates 
ak,bk and (in the periodic case) on the additional parameter A: 

N N N 

T{a, b,X) = J2 hEkk + A X! Ek+i,k + J2 o-kEk,k+i- (2.7) 

k=l k=l k=l 

Here Ejk stands for the matrix whose only nonzero entry on the intersection of the j'th row 
and the kth column is equal to 1. In the periodic case we have Ejy+i^N = -E'i.at, Ejy^^^i = 
Ejy^, in the open-end case we set A = 1, and i^Ar+i^Tv = -E'Af.iv+i = 0. 
The flow ( ^.11 ) is equivalent to the matrix differential equation 

f=[T,B], (2.8) 

where 

N N 

Bia,b,X) = J2bkEkk + \Y.Ek+i,k, (2.9) 

k=l k=l 



and the map (2]^) is equivalent to the matrix difference equation 

f = B^TB, (2.10) 
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where 

N N 

B(a, b,\) = Y^ PkEkk + hXY^ Ek+i,k. (2.11) 

k=l k=l 

The spectral invariants of the matrix T(a, 6, A) serve as integrals of motion for the flow 
(^■11), as well as for the map (|2.4|). In particular, 



An involutivity of the spectral invariants of the matrix T follows from the r-matrix 
theory. Namely, it turns out that the matrices T(a, 6, A) form a Poisson submanifold for a 
certain linear r-matrix bracket on the algebra of N x N matrices (in the open-end case) 
or on the loop algebra over the previous one (in the periodic case) (cf. [§), and (|2.2|) 
gives exactly the corrdinate representation of this r-matrix bracket restricted to the set 
of the matrices T(a, b, A). 

3 Continuous and discrete time Toda flows 
in a new parametrization 

The usual parametrization of the variables (0^,6^.) by the canonically conjugated ones 
{xk,Pk) is: 

Qk = exp(a;fc+i - Xk), bk = Pk- 

This corresponds to the bracket ( |2.2| ) and results in the most usual form of the Toda 
lattice and in its discretization studied in [0. 

Here we chose a different parametrization, which, however, corresponds to the same 
bracket ( p.2|) and is in a sense dual to the parametrization just mentioned: 



ttk = exp(pfc), bk = Xk- Xk-i. (3.1) 
The corresponding Hamiltonian function takes the form 

N N 

Ht{x,p) = J2 exp(pfe) + - ^(a;fc - Xk-if, (3.2) 

k=l ^ fc=l 

which generates the canonical equations of motion 

Xk = dHT/dpk = exp{pk), 

Pk = -dH-r/dxk = Xk+i - 2xk + Xk-i- 
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These equations of motion imply the Newtonian ones ( |1.2|) . 

We turn now to a less straightforward case of discrete equations of motion. 

Theorem 1. In the parametrization ( p.l|) the equations of motion (|2.4| ) may be pre- 
sented in the form of the following two equations: 

h exp(pfc) = (xk - Xk) (l + h{xk - , (3.3) 

h exp(pfc) = {xk - Xk) (l + h{xk+i - 2^)) , (3.4) 



which imply also the Newtonian equations of motion (|L 

Proof. The second equation of motion in (|2.4| ) together with the second equation in 
( p.l|) implies that 

hak ^ hok-i 



Xk Xk „ Xk—l Xk—l 



Pk Pk-1 

Hence the expression on the left-hand side does not depend on k, and setting it equal to 
0, we get: 

hok ^ , . 

Xk-Xk- (3.5) 



Substituting this into the recurrent relation (p.5|), we get: 



= 1 + /i(a;fc - Xfe_i). (3.6) 



The last two formulas imply the expression ( |3.3|) for hok = hexp{pk). Finally, the first 
equation of motion in (|2.4|) implies the expression ( p.4|) , if one takes into account 



4 Simplest flow of the Bruschi— Ragnisco hierarchy 

This section is devoted to reminding the known facts about the simplest flow of the 
Bruschi-Ragnisco hierarchy 0,0. 

It is given by the equations of motion: 

bk = bk+iCk - bkCk-i, Ck = Ck{ck - Ck-i). (4.1) 

It may be considered either under open-end boundary conditions (cq = bjy+i = 0), or 
under periodic ones (all the subscripts are taken (mod A^), so that Cq = cn, b^+i = &i)- 
The flow ( [4. ID is Hamiltonian with respect to the following Poisson bracket: 

{ck, bk} = -{ck, bk+i} = Ck, (4.2) 
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with the Hamiltonian function 

N 

Hbr = h+iCk- (4.3) 

k=l 

The Lax representation of the flow ( |4.1|) is given in terms of the N x N rank 1 Lax 
matrix T{b, c) with the entries given by: 

' bj UlZl c, , k<j 
Tkj = ^ (4.4) 

. bj c,) , k> J 

The flow ( ^.1| ) is equivalent to the following matrix differential equation: 

f=[T,M], (4.5) 

where M is a constant matrix, so that the Lax equation is in fact linear. For example, in 
the open-end case 

N-l 

M=Y Ek+i,k- (4.6) 

k=l 

Hence the solution of the flow ( [4.1|) can be read off the following explicit formula: 

T{t) = exp(-tM)T(0) exp(tM). (4.7) 

The integrals of motion for the flow are given not by the spectral invariants of 

the matrix T{b,c) (they are all dependend because this matrix has reank 1), but by the 
following set of linear functionals on the Lax matrix: tr (M'^T). In particular, 

HBK = tr{MT). 

The involutivity of such linear functionals on the matrix T with respect to the bracket 
( [4.2|) follows from its interpretation as a coordinate representation of a standard Lie- 
Poisson bracket on the algebra gl{N), restricted to a Poisson submanifold of rank 1 ma- 
trices [Q. Taking such functions as Hamiltonians, i.e. replacing in the previous formulas 
M by M™, l<m<A^ — 1, we get the whole Bruschi-Ragnisco hierarchy. 

According to the philosphy of 0, 0, the discrete time Bruschi-Ragnisco lattice be- 
longs to the same hierarchy, i.e. it shares the Lax matrix with the continuous time one, 
and its explicit solution should be given by 

T{nh) = {I + /iM)^"T(0)(/ + /iM)". (4.8) 

Hence the corresponding discrete Lax equation should have the form 

f = {I + hM)-^T{I + hM). (4.9) 
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A direct calculation shows that the last equation in the coordinates (6^, c^) has the 
form: _ 

bkil + hck-i) = hk + hbk+iCk, 4 = Ck -—rz — ■ (4.10) 

1 + hck-i 

By construction, this map is Poisson with respect to the bracket ( [4.2|) . 

5 Continuous and discrete time Bruschi— Ragnisco 
lattices in the Newtonian form 

Observe now that the bracket ( [4.2|) for the Bruschi-Ragnisco lattice formally is identical 
with that for the Toda lattice (|2.2| ). 
This suggest the parametrization 

Ck = exp(pfc), bk = Xk- Xk-i- (5.1) 

The corresponding Hamiltonian function H^r takes the form 

N 

Hbr{x,p) = ey:p{pk){xk+i - Xk), (5.2) 
fc=i 

which generates the canonical equations of motion 

Xk = dHBn/dpk = ey:p{pk)ixk+i - Xk), 
Pk = -dHBR/dxk = exp(pfc) - exp{pk-i). 

These equations of motion imply the Newtonian ones ( |1.3D . 

Turning to a less straightforward case of discrete equations of motion, we have: 
Theorem 2. In the parametrization (|5.1| ) the equations of motion (|4.1CI| ) may be 

presented in the form of the following two equations: 

, / N {Xk- Xk) {Xk - Xk-l) 

/lexp(pfe) = -, (5.3) 

[Xk - Xk-l) [Xk+l - Xk) 

heMPk) = 7^^^^^, (5.4) 

[Xk+l — Xk) 

which imply also the Newtonian equations of motion (|1.6D. 

Proof. We start with the chain of identities following from the equations of motion 
( [4.1(j| ), independently of any parametrization: 

^ bk+i _ hbk+iCk _ hbk+iCk _ ^ bk 

''^^ 1 + hck 1 + hck 1 + hck-i ^ 1 + hck-i 
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Now we substitute 6^+1 = Xk+i — into the leftmost part of this chain, bk = Xk — x^-i 
into the rightmost part, and arrive at the conclusion that 



1 + hck 1 + /iCfc-i 

So, the expression on the left hand side of this equality does not depend on k. Setting it 
equal to 0, we obtain: 

, , , ~ h+i Xk+i - Xk 
l + hck = — = —, 

Xk+l ~ Xk Xk+1 — Xk 

which is, clearly, equivalent to ( ^.41) , and, together with the second equation of motion in 
( CT ), implies also (p]). ■ 

6 Simplest flow of the relativistic Toda hierarchy 

In this section we recall some known facts about the relativistic Toda lattice. The sources 
are: |^ and the references therein. 

The simplest flow of the relativistic Toda hierarchy is: 

dk = dk{ck - Ck-i) , Ck = Ck{dk+i + Ck+i - dk - Ck-i) . (6.1) 

It may be considered either under open-end boundary conditions {d]\f+i = Cq = cn = 0) , 
or under periodic ones (all the subscripts are taken (mod A^), so that dj\f+i = di, Cq = cat, 

Cn+1 = Ci). 

It is Hamiltonian with respect to the Poisson bracket 

{cfc,4+i} = -Cfc, {ck,dk} = Ck, {4,4+i} = Cfc, (6.2) 
with a Hamiltonian function 

^ N N 

Hrt = ^ XK'^fc + Ck-if + '^{dk + Ck-i)ck. (6.3) 

^ fc=l k=l 

An integrable discretization of the flow ( |6.1j ) derived in is given by the difference 
equations 

J J Ofc+i ~ hdk+i _ afc+i + hck+i ,„ 

dk = dk — , Ck = Ck — , (6.4) 

ttfe — ndk ak + iiCk 

where = afc(c, d) is defined as a unique set of functions satisfying the recurrent relation 

ttfc+i = 1 + hdk+i + — (6.5) 



together with an asymptotic relation 

afc = l + /i(4 + Cfe_i)+0(/i2). (6.6) 



In the open-end case, due to Cq = 0, one obtains from ( |6.5|) the following finite continued 
fractions expressions for a^: 

hci 

ai = 1 + hdi] a2 = 1 + hd2 ' 



= 1 + hdi\f + 



1 + hdi' 
hcN-i 



1 + hdN-i H 



1 + hdN-2 + 



hci 

+ 



1 + hdi 

In the periodic case ( |6.5|) , (|6.6| ) uniquely define a^'s as A^-periodic infinite continued 
fractions. It can be proved that for h small enough these continued fractions converge 
and their values satisfy ( |6.6|) . 

The map (|6.4|) is Poisson with respect to the bracket (|6.2|). 



Both the continuous time flow (|6.1|) and the discrete time one (|6.4|) admit Lax repre- 
sentations with the same Lax matrix, given by 

T(c, d, A) = L(c, d, X)U-\c, d, A), (6.7) 

where L and U are two N x N matrices depending on the phase space coordinates c^, dk 
and (in the periodic case) on the additional parameter A: 

N N 

L{c,d,X) = Y.dkEkk + XT.Ek+i,k, (6.8) 

k=l k=l 

N N 

U{c,d,\) = Y.Ekk- \~'T.^kEk,k+i. (6.9) 

k=l k=l 

The flow ( |6.1D is equivalent to the matrix differential equation 

f=[T,A], (6.10) 

where 

N N 

A{c, d, A) = ^(4 + Ck-i)Ekk + A ^ Ek+i,k. (6.11) 

k=l k=l 
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The map ( |6.4|) is equivalent to the matrix difference equation 

T = A-^TA, (6.12) 

where 

N N 

A(c, d,X) = Y^ akEkk + h\Y^ Ek+i,k, (6.13) 

fc=l k=l 

and the quantities are defined by the recurrent relations ( |6.5|) . 

The spectral invariants of the matrices T(c, d, A) serve as integrals of motion for the 
flow ( |6.1| ), as well as for the map ( p.4| ). In particular, 

HnT = ItiiT'). 

As in the case of the Toda lattice, the involutivity of the spectral invariants of the 
matrix T with respect to the bracket { p.'2\) follows from the r-matrix theory. Namely, 
it turns out that the matrices T(c, d, A) form another Poisson submanifold for the same 
bracket on the matrix algebra as in the Toda case, and (|6.2|) serves as a coordinate 
representation of this bracket restricted to this manifold (cf. [0). 



7 Continuous and discrete time relativistic Toda 
lattices in a new parametrization 

It is easy to see that the following parametrization of the variables {ck, d^) by the canon- 
ically conjugated variables {xk,Pk) leads to the linear bracket 



dk = Xk- Xk-i - exp(pfc-i), Ck = exp{pk). (7.1) 

Let us find out how do the equations of motion look in this parametrization. We start 
with (13). 

Obviously, the function if^x takes the form 

N ^ N 

Hrt = exp(pfc)(a;fc - Xk-i) + Xk-i^- (7.2) 

k=l ^ k=\ 



Correspondingly, the flow (6J-) takes the form of canonical equations of motion: 
Xk = dHBT/dpk = exp{pk){xk-Xk-i), 

Pk = -dHKr/dxk = Xk+i - 2xk + Xk-i + expipk+i) - exp{pk). 
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As an immediate consequence of these equations one gets the Newtonian equations of 
motion ( |1 . 1| ) . 

We turn now to a less straightforward case of discrete equations of motion. 



Theorem 3. In the parametrization (TJ.) the equations of motion (|6.4| ) may be pre- 
sented in the form of the following two equations: 

hexp{pk) = , f''^ . (l + h{xk - Xfc-i)), (7.3) 

[Xk - Xk-l) ^ ^ 

hexpiPk) = 7- ^ 7^ ^ r (1 + h{xk+i -Xk)). (7.4) 



which imply also the Newtonian equations of motion ( |1.4]) . 

Proof. We start with derivation of several useful formulas which do not depend on 
the parametrization of the variables {ck,dk). 

The equations of motion ( |6.4|) together with the recurrent relation ( |6.5|) imply: 



7 , ~ / , , ifc A - hdk+i , , 

dk + Ck-i = \ dk + Ck-i 7-1 • (7-5) 



But, applying once more (p. 5]), we derive 



hdk + hck-i = hdk + ak{ak-l - hdu) = (a^ - l)(afc - hdk), 

SO that ( |7.5| ) may be rewritten as 

h{dk + Cfc-i) = {ak - l)(afc+i - hdk+i). (7.6) 
It is easy to see that, according to ( |675D , 

(cifc — l)(aA:+i — hdk+i) = h{dk+i + Ck) — (ofc+i — afc) , 
so that ( |7.6| ) implies: 

afc+i - ak = h{dk+i + Ck) - h{dk + Ck-i). (7.7) 

Now we are in a position to make use of the concrete parametrization (|7.1|) . Namely, 
due to dk+i + Cfc = Xk+i — Xk the last formula implies that 

afc - h{xk - Xk-i) = ak+i - h{xk+i - Xk), 

i.e. the expression on the left-hand side does not depend on k. Setting this equal to 1, 
we arrive at: 

ak = l + h{xk-Xk-i). (7.8) 
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Now the formula (|6.5| ) upon use of (|7.1| ), { \l.8\) may be presented as 

1 + h{xk+i - Xfc) = 1 + h[Xk+i -Xk- exp{pk) ) + ^ , , . r, 

V ^ l + h{xk-Xk-i) 

which can be easily resolved for exp{pk) and gives (iTiSl) . Next, O together with the 
expression ( |7.3| ) for = exp(pfc) implies 

ak + hCk = 7 r( 1 + h{xk - Xk-i)). 

{Xk - Xk-l) ^ ' 

Substituting this into the second equation of motion in (p.4[), we finally obtain ( |7.4| ). ■ 

8 Conclusion 

The relations between three well known lattice systems, established in the present paper, 
may be symbolically presented as 

RTL = TL + BRL , dRTL = dTL x dBRL . 

They are, in my opinion, rather beautiful and unexpected. However, they probably are 
nothing more than a pure curiosity, which demonstrates once more: the field of integrable 
systems of classical mechanics, even in its best studied parts, is far from being exhausted. 
Many important and less important but funny findings await us in this fascinating world. 
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